Mon. Not. R. Astron. Soc. 000, 000-000 (1995) 



Recovering the Intrinsic Metallicity Distribution of 
Elliptical Galaxies 

L. Ciotti 1,2 , M. Stiavelli 1 ' 3 and A. Braccesi 4 

1 Dept. of Theoretical Physics, Oxford University, UK 

2 Osservatorio Astronomico di Bologna, 140100 Bologna, Italy 

3 Scuola Normale Superiore, 156126 Pisa, Italy 

4 Dipartimento di Astronomia, Universita di Bologna, 140100 Bologna, Italy 



22 May 1995 



ABSTRACT 

We address the problem of deriving, from the observed projected metallicity gradients, 
the intrinsic metallicity distribution of elliptical galaxies as a function of their integrals 
of motion. The method is illustrated by an application to anisotropic spherical Hernquist 
models. We also compare the derived metallicity distribution with those expected from 
two very simple models of galaxy formation and find that the more dissipative scheme 
agrees better with the typical metallicity distribution of ellipticals. 
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The aim of this paper is to give a simple recipe to derive the 
metallicity distribution of galaxies as a function of their inte- 
grals of the motion. Elliptical galaxies are known to possess 
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"■"^metallicity gradients (Peletier 1989, Gorgas, Efstathiou and 
r^H Salamanca 1990, Davies, Sadler and Peletier 1993, Carollo, 
^■Danziger and Buson 1993), established mainly by measuring 
Q the Mg2 line strength, which is one of the most prominent 
?H features in the optical spectra of early type galaxies and K 
^3 and M type giants. The observed metallicities, as determined 
by this line frequently show variations of a factor 2 from the 
center to the effective radius. The observed gradients are the 
result of projection and orbital mixing. The former is easy 
to model and was in fact shown by Peletier (1989) to lead 
to only minor changes in the observed metallicity gradients. 
The latter is more complex and has not been studied in de- 
tail. Orbital mixing arises because stars at a given radius 
may have their apocenters spanning a wide range of radii. 
Thus, the mean metallicity at any given point inside the 
galaxy is the result of the metallicity distribution in phase 
space weighted by the galaxian distribution function. 

In order to investigate these effects, we develop a tech- 
nique to derive the metallicity distribution in phase space 
from the observed metallicity gradients of for spherical 
galaxies. The technique is a simple extension of Eddington's 
inversion to derive the distribution function from the galaxy 
density profile. The method is described in Sect. 2 and 3, and 
applied to anisotropic spherical Hernquist models in Sect. 4. 
In Sect. 5 we compare the derived phase space metallicity 
distribution to the one predicted by simple models of galaxy 



formation (Sect. 5), finding the simple prediction of dissipa- 
tive models in better agreement with observations. 



2 FROM Z P TO Z D 

For a spherical galaxy the luminosity density v(r) is ob- 
tained by deprojecting its surface brightness profile I(R). 



v(r) 



dl(R) dR 



dR V-R 2 - r 2 



(1) 



where r is the spatial radius, and R the projected one (see, 
e.g., Binney and Tremaine 1987, hereafter BT). The ob- 
served (projected) metallicity Zp(R) is the result of a projec- 
tion along the line-of-sight of the luminosity- weighted spatial 
mean metallicity Zr>{r): 



I(R)Z P (R) = 2 



Jr. 



j(t)Zd (r)rdr 

which can be inverted to yield: 

dI(R)Z P (R) 



v(t)Zb(t) 



dR 



dR 



V^R 2 " 



(2) 



(3) 



Note that Zn(r) as derived from the ratio of Eqs. (3) and 
(1), does not depend on the unknown galaxy mass-to-light 
ratio T(r) = p(r)/v(r). 



3 RECOVERING Z(Q) 

Assuming that the previous galaxy is described by a distri- 
bution function / = f(Q), the relation between its density 
p and / is 
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pq{ t ) = p( r ) 1 + 



47T 



/ f(Q)V' 

Jo 



2(*-Q)dQ, (4) 



where Q = £-L 2 /2r 2 a , £ = <£(r)-v 2 /2 is the stellar binding 
energy per unit mass, ^(r) is the gravitational potential, I? 
is the square modulus of the angular momentum, and r a 
is the so-called anisotropy radius (Osipkov, 1979; Merritt, 
1985; see also BT for a general account). For r a = oo the 
distribution function depends on the energy only, and the 
velocity dispersion tensor is globally isotropic. For r a = all 
the stellar orbits are radial. Equation (4) can be inverted, 
namely (see BT): 



f(Q) 



VStt 2 dQ J Q <2* ^/^^ ^ *" 



(5) 



Therefore, once ~£(r) is known, one can formally derive, 
for any fixed anisotropy radius, an equilibrium distribution 
function f(Q) for a spherically symmetric system with a 
given 1(H)- 

The quantity Zn(r) is the mean metallicity of all stars 
observed at the radius r; these stars can be characterized by 
very different Q values, i.e. their orbits lie in different regions 
of the phase space. In order to derive the intrinsic metallic- 
ity distribution in phase space, the function Z(Q), one pro- 
cedes in a similar way that in Eqs. (l)-(3). The metallicity 
dependence of the mass-to-light ratio of a stellar population, 
T = ¥(Z), gives the two equations: 
/»* 

vq{t) = ±-k f{Q)T-\Z)^2{¥~Q)dQ, (6) 
Jo 

and 

/»* 

v Q {r)Z D {r) = 4tt / f(Q)T-\Z)Z{Q)y/2{* -Q)dQ. (7) 

Jo 

Their functional inversion is the same as in Eq. (5): 
1 d f Q dug 



f{Q)1-\Z) 



I 



VsV dQ J VQ- *' 



and 



f{Q)T-\Z)Z{Q) 



V^tt 2 dQ 
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dvqZn d*S! 



d<5> V^T 7 *' 



(8) 



(9) 



The functions needed for the inversion are provided by Eq. 
(1) and Eq. (3), and we can derive separately the product 
f(Q)T~ 1 (Z) and the intrinsic metallicity Z{Q), which does 
not depend on the knowledge of T. Moreover, if we can de- 
termine T from synthetic stellar population, e.g., for a 12 
Gyrs population: 



log 10 Ty ~ 1.178 + 0.1927 log 10 Z 



(10) 



(Peletier, 1989), then from Eq. (8) we can obtain also /, 
and from that the projected velocity dispersion ap(R). A 
difficulty may arise if the derived Z(Q) and f(Q) are nega- 
tive for some permitted value of Q. Indeed, it is well known 
that the previous inversion procedure does not guarantee 
the positivity of the recovered function. 



4 A TOY MODEL 

We have applied the method to the case of Hernquist mod- 
els (Hernquist, 1990), that approximate the de Vaucouleurs 



Figure 1. Left panel: projected (solid line) and deprojected (dot- 
ted line) metallicity profiles as a function of radius. Right panel: 
the derived intrinsic metallicity profile is shown, for the isotropic 
model (solid line), and for anisotropic model with r a = R e . The 
two short-dashed curves refer to the dissipationless galaxy forma- 
tion models described in Sect. 5, and the straight, long-dashed 
line refers to a simple dissipative formation model. 



il 1 / 4 law, and for which many properties can be evaluated 
analytically. The model investigated is characterized by a 
total mass 10 11 Mq, an effective radius R e = 4 kpc, and a 
projected metallicity profile given by 



Z P (R) 



Zp(0) 



0.18 1og 10 (l + jjj) 



,(H) 



which agrees with the observations of Mg2 gradients in el- 
lipticals calibrated to metallicity by using equation (1) of 
Bender, Burstein and Faber (1993) for a population 12 Gyrs 
old. Although steeper gradients will have a larger effect, we 
consider this slope to be representative of the mean gra- 
dients measured in ellipticals. Inside Rz = -R e /30, Zp(R) 
is taken to assume the central value Zp(Q) = 2Zq, where 
Zq ~ 0.02 is the solar metallicity. Our results are shown in 
Fig. 1, where in the left panel we show the projected and 
the deprojected metallicity profiles, Z/Zq, as a function of 
radius. In the right panel we show the intrinsic metallic- 
ity of stars as a function of Q, for a globally isotropic model 
and for an anisotropic model characterized by r a = Re - Note 
that the phase-space metallicity distribution obtained is pos- 
itive everywhere, and, for the isotropic case, nicely close to 
a straight line over a wide interval in Q. The same linear 
behavior is obtained also for steeper and flatter metallicity 
gradients (not shown in Fig. (1)). In the next section we will 
discuss some implications of this result. 



5 COMPARISON WITH FORMATION 
MODELS 

The formalism derived above also allows one to set some 
constraints on galaxy formation models. Rather than con- 
sidering a given galaxy, here we focus on the model described 
above and taken to be representative o{ typical ellipticals. As 
a first comparison of the derived intrinsic metallicity Z(£) 
(here we consider only the isotropic case, where Q reduces 
to the binding energy £) with formation models, we have 
considered two different, extreme but simple, scenarios. The 
first scenario is based on the assumption that galaxies col- 
lapse dissipatively and, after they have essentially settled 
on their final density profile, form stars. The binding en- 
ergy distribution in such a system would match closely the 
gravitational potential. In this scenario one would expect 
Z(£) cc £. This is the underlying assumption leading several 
authors to compare metallicities with escape velocities (see, 
e.g., Franx and Illingworth 1988, Davies et al. 1993, Carollo 
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and Danziger 1994). Despite the simplicity of the argument 
and the lack of detailed modelling, these authors confirmed 
that a good correlation is indeed observed between (pro- 
jected) metallicities and escape velocities averaged along the 
line of sight. 

The second scenario is based on the assumption that 
galaxies collapses without dissipation after having formed 
their stars on timescales shorter to dynamical ones. We as- 
sume that stars form out of a top hat gas cloud of constant 
density po . Such a model has a gravitational potential given 
by: 



*(0 = *(0)(1 - k 2 ), 



(12) 



with Vt(0) and £ being, respectively, the central value of the 
gravitational potential and the radius relative to the outer 
boundary of the top hat disturbance. 

At each radius, star formation and enrichment proceed 
following the simple model illustrated in, e.g., BT, until the 
gas becomes unbound due to SN energy release. We assume 
that each radial shell evolves as a closed box model. The 
metallicity profile of stars in the (pre-collapse) protogalaxy 
is then given, in terms of the gas metallicity, by: 



z*{£) = p 



1 - [1 + Z gas ($)/p] e xp[-Z gas ($)/p] 

1 - exp -[Zgas{t)/P] 



(13) 



where p is the chemical yield. In the limit of very large 
Zgas/p Eq. (13) gives = p. In the opposite limit, 

Zgas/p <. 1 one has Z* = Z gas /2, which implies p > 2Z*. 
Since the element production of the supernovae is propor- 
tional to their energy release, and, at the onset of the galac- 
tic wind, the energy release per unit mass equals the local 
gravitational potential, one has: 



Zgaa{i) — Z gas (0) 



mi 

*(0)' 



(14) 



Note that the stars that have formed do not acquire 
significant velocities while star formation is active, i.e., 
£ = VP (£), since the basic assumption of the dissipation- 
less collapse scheme is that the star formation time scale is 
shorter than the dynamical time scale. Assuming that the 
evolution of gravitational potential during the wind phase 
can be neglected, the stellar density, after the wind has oc- 
curred, is obtained from Eq. (9-21) of BT, namely: 



p*(0 

Po 



1 — exp 



Z 9 as(0) *(Q' 



P 



*(0) 



(15) 



where Z gas (0) is the central gas metallicity. This allows us 
to compute the cumulative energy distribution at the occur- 
rence of the wind, i.e. the mass fraction of stars more bound 
than a given energy. By assuming that the subsequent dis- 
sipationless collapse does not alter this distribution, as ac- 
tually verified in N-body simulations, one predicts the final 
Z(£). This is accomplished by remapping the initial cumula- 
tive energy distribution into that of the an Hernquist sphere 
with the same total mass. 

The predictions from these two models are shown in Fig. 
1 (right panel, dashed lines) where they are compared to the 
intrinsic metallicity derived by our inversion procedure. The 
two short dashed lines represent the top hat gas cloud model 
prediction for the intrinsic metallicity of a galaxy with a final 
central star metallicity of 2Zq. The two lines correspond to 



the minimum value of p giving the desired central stellar 
metallicity, p = 0.04, and to p = 0.1. All curves for p 3> 0.04 
are similar. Note that the simple dissipative scheme (Fig. 
1, long dashed line) appears to better support the typical 
observed gradient than the simplified dissipationless scheme. 



6 CONCLUSIONS 

We have proposed a simple inversion procedure allowing one 
to derive the dependence of metallicity on the integrals of 
motion for a spherical galaxy, under the assumption of a 
Osipkov-Merritt orbital anisotropy distribution. Although 
applied only to an analytical Hernquist model, the technique 
described here could be applied to real (close to spherical) 
galaxies for which surface brightness and metallicity pro- 
files are available. The ideas of this paper can be general- 
ized to oblate two-integrals models constructed following the 
Hunter fc Qian (1993) technique. The intrinsic metallicity as 
a function of binding energy and angular momentum would 
then be obtained by inversion from a two-dimensional Mg2 
map. 

We have also attempted a cursory discussion of the pre- 
dictions of two simple galaxy formation scenarios for the 
intrinsic metallicity distribution as a function of the bind- 
ing energy. We conclude that for typical ellipticals the more 
dissipative scheme is to be preferred to the dissipationless 
scheme. In order to derive more firm conclusions, the metal- 
licity derived from gas dynamical simulations should be con- 
sidered. 
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